Bias-induced insulator-metal transition in organic electronics 
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We investigate the bias-induced insulator-metal transition in organic electronics devices, on the 
basis of the Su-SchriefFer-Heeger model [W. P. Su, J. R. Schrieflter, and A. J. Heeger, Phys. Rev. 
B 22, 2099 (1980)] combined with the non-equilibrium Green's function formalism. The insulator- 
metal transition is explained with the energy levels crossover that eliminates the Peierls phase 
[R. Peierls, Quantum Theory of Solids, Oxford University Press, Oxford, 1955] and delocalizes the 
electron states near the threshold voltage. This may account for the experimental observations on 
the devices that exhibit intrinsic bistable conductance switching with large on-off ratio. 

PACS numbers: 71.30.-|-h, 72.80.Le, 73.63.-b 
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In recent years molecular electronics, which uses in- 
dividual molecules as possible electrical switches, has 
grown up rapidly since the discoveries of negative dif- 
ferential resistance and bistable conductance switching 
(CS) [J, [I]- Individual molecules can be wired to serve 
as diodes, transistors, and other electronic devices at the 
heart of computer chips [1, 0) S • For the next gener- 
ation of technology, molecular electronics should exhibit 
a large on-off ratio (>50:1). It should also be intrinsic 
and controllable, rather than artificial due to for example 
metal filament [7.] . Some mechanisms such as the confor- 
mational change 0, charging effect and polaron mod- 
els [1, II] have been proposed. However, none of them 
can fully account for the large intrinsic on-off ratio of 
molecular switching observed in experiments [lol . [ill . . 
In this work, we shall propose another possible mecha- 
nism, in which the intrinsic CS with large on-off ratio is 
controlled by the bias-triggered Peierls phase elimination 
and energy levels crossover. 

Let us start with the following observations. Most 
molecular devices use short-chain 7r-conjugated organic 
molecules as the active region, and are fabricated as 
metal-molecule-metal sandwiched structures. A chain- 
like nanoscale organic molecule belongs to the fam- 
ily of organic semiconductors (OSEs). Physically, an 
OSE is charactered by the Peierls instability; i.e., 
the strong electron-phonon (e-ph) interaction that in- 
duces the single-double bond alternation and a relative 
large energy gap between highest-occupied and lowest- 
unoccupied molecular orbital (HOMO and LUMO) [11]. 
It is the Peierls instability that largely localizes the 
TT-electrons of carbon and contributes to the low con- 
ductance of OSEs, no matter they are fabricated via 
the bottom-up self-assembled monolayer or the top- 
down lithography technique. Apparently, delocalizing tt- 
electrons will improve the conductivity of OSEs. One 
method is the creation of nonlinear excitations (solitons 
or polarons) to act as charge carriers [l3]- Another ap- 
proach is to eliminate the Peierls phase directly. Increas- 
ing temperature may melt the Peie rls p hase, leading to 
an insulator-metal (I-M) transition |13j . 



An effective and controllable method in electronics 
should however be the one triggered by the bias or 
other external control field. In this connection and 
also in contact with some typical experimental systems 
(Til , we present in this work a theoretical study on 
the bias-induced I-M transition in a model nanoscale 
metal/OSE/metal sandwich structure. We will elucidate 
the fact that the bias-triggered energy levels crossover 
and its resulting Peierls phase elimination can largely ac- 
count for the observed CS. 

We combine the celebrated Su-Schrieffer-Heeger (SSH) 
model ]14] and the nonequilibrium Green's function 
(NEGF) formahsm [H, [3 to study organic electronics. 
For the spin-independent charge transport considered in 
this work, the SSH Hamiltonian for the OSE electrons 
coupled adiabatically with the lattice displacements as- 
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Here, (c„) denotes the creation (annihilation) opera- 
tor of an electron at the n*'' site of the OSE, and the 
on-site energy, to the zero-displacement hopping integral, 
and ti the nondegeneracy parameter, respectively. The 
lattice distortion is written in terms of the bond distances 
{Un — Wn+i —Un} and treated classically, with the spring 
constant Ko and the adiabatic e-ph coupling constant ao- 
The L and R electrodes are chosen to be an identical non- 
magnetic metal. They are individually described by the 
one-dimensional single-band tight-binding Hamiltonian, 
being of the on-site energy e / and the nearest neighbor 
transfer integral t f . 

At the mean-field level of correlations, one can evaluate 
the current based on the NEGF formalism [Tl, [3 , 
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Here, f{E, /^l/r) is the Fermi distribution function at the 
lead chemical potential /^l/r- The trace term in the inte- 
grand is the transmission coefBcient function T{E). The 
involving Ti^/^ denotes the lead reservoir-induced broad- 
ening matrix, while C and — (G" )^ are the retarded 
and the advanced single-particle Green's functions for the 
central scattering region (S-region), under a finite applied 
bias. This region consists of the OSE together with a 
number of metal atoms attached to each of its ends. For 
a quantum open system, the effects of both the charging 
from the electrodes and the external potential from the 
bias voltage are contained in the nonequilibrium reduced 
density matrix, p = ^ J2n=h.K J G'TaG'^f{E,^a)dE, 
of the coupled e-ph system [15{ . The reduced density 
matrix should be evaluated in a self-consistent manner, 
together with the Poisson's equation and the Hellman- 
Feynman (H-F) variation theorem for lattice distortion: 
9[Tr(iJsp)]/9M„ — 0. The self-consistent procedure as 
well as the time-to-voltage mapping method to handle 
the time-dependent bias sweeping have been described 
in detail elsewhere [1, [l3| • 

In the following calculation of the hysteretic I — V 
characteristics, we choose a constant rate of bias sweep- 
ing to be 0.1 V/sec, as both the numerical efficiency 
and the typical experimental value being of 0.01 ~ 0.1 
V/sec are concerned. Note also that SSH-I-NEGF formu- 
lation includes two main approximations. The first one 
is the static lattice description involving in the model 
Hamiltonian and the aforementioned H-F variation the- 
orem. This approximation is valid when the lattice fluc- 
tuation time scale is longer than the inverse coherent 
charge transfer coupling. The second approximation is 
the meanfield treatment of electron-electron correlation 
that involves in the NEGF formulation and also the 
model Hamiltonian. This is a common approach, espe- 
cially in treating strong e-ph systems. 

We choose the SSH parameters in Eq. ([T]) the typ- 
ical values: to = 2.5 eV, h = 0.04 eV, ao = 4.2 
eV/A, and Ko = 21.0 eV/A^. The transfer integral for 
metal electrodes is tj = 3.0 eV to produce a wide band 
i^tf), and the metal-OSE coupling transfer integral is 
^metal — OSE — 0.5{to + tf). Without losing generality, we 
set Ep — eo — £f — 0, where Ep is the equilibrium Fermi 
level at the zero bias V = 0. 

We then calculate the nonlinear current [Eq. ^] as a 
function of sweeping bias of the model metal/OSE/metal 
electronics at T = 11 K by using the time-to- voltage 
mapping method [l7j. The resulting I-V curve, pre- 
sented in FiglU clearly shows hysteretic CS characters. 
Note that the model OSE is of No = 40, which is about 
5nm and comparable to the experiments in Refs. [13 and 
[Tsl . More than 30 years ago, Simmons et al. reported 
an electroformed metal-insulator-metal diode with re- 
versible voltage memory effects [l9|. Recently, Bozano 
et al. had also demonstrated a similar resistive switch- 
ing phenomenon in OSE layers [20| . They had proposed 
several necessary features of an operational switch, as 
summarized in Ref. '20. In comparison with their exper- 



iments, especially the Fig. 3 in Ref. [l^, our model calcu- 
lation dose possess most of the necessary features they 
proposed. These include (1) a distinct bistable on/off 
state corresponding to low/high resistance with an on- 
off ratio > 20:1; (2) a local maximum in current /max 
which nonlinearly depends on the size of OSE and the 
metal-OSE coupling (varying from several nA fiA) at 
voltage VJnax (~ 4.5 V); (3) a threshold of voltage Vth (~ 
3.7 V) at which the high-current (switch-on) state is es- 
tablished; (4) reproducible switching and reading which 
can be achieved with applied positive or negative volt- 
age. The insert in FigU] shows the nonlinear dependence 
of the threshold voltage Vth on the OSE chain length 
No- This confirms the generality of our results in organic 
electronics devices with different sizes. 

Let us discuss the mechanism of the bistable CS, as in- 
dicated in FiglU The OSE substructure is found charge 
neutral, indicating there are no charging and polaron 
effects below the saturation voltage. The abrupt in- 
crease of current near the threshold Vth at 3.7 V should 
therefore be interpreted with an alternative mechanism. 
First of all, we shall show that the observed abrupt cur- 
rent change relates to an I-M transition. Presented in 
Figl2] are the changes of transmission coefficient T{E) 
and lattice distortion ?/„ around the threshold voltage. 
Evidently, the Peierls phase of the OSE substructure is 
largely eliminated near Vth- The key features in Fig[2] 
include (1) the rapid increasing of transmission coeffi- 
cient in the window between the indicated /il and /iR in 
FigEl^a); (2) the distinct reducing of the lattice distortion 
or single-double bond alternation [Figl^l^b)]; and also (3) 
the vanishing of the energy gap and the merging of va- 
lence and conductance bands, which will be explained 
later, cf. Figl3]Ja). The above features all indicate an I- 
M transition around Vth. When the bias sweeps down, 
the Peierls phase, as indicated by the lattice distortion in 
FigHfb) (see also FigO for energy gap) , is reformed re- 
sulting in a metal-insulator (M-I) transition around 2.8 
V. The delayed M-I transition, in comparing with the 
I-M transition voltage, is due to the lattice distortion re- 
organization, resulting in the hysteretic loop (~ 1 V) in 
the I-V curve in FiglTJ 

To further elucidate the bias-induced I-M/M-I transi- 
tion mechanism, we also calculate the 'reduced' energy 
levels and wave functions of the OSE substructure by di- 
agonalizing its Hamiltonian. It consists of Eq. ([T]), with 
the nonequilibrium lattice parameters being evaluated 
in the self-consistent manner as described earlier. The 
Hamiltonian contains also the electric potential variation 
in the OSE structure, which had also been determined 
self-consistently. The informations on this Hamiltonian 
are shown to be able to reflect the bias-induced intrinsic 
changes of the OSE in the device. 

The resulting energy levels from the aforementioned 
self-consistent Hamiltonian, as they vary with the up- 
sweeping and down-sweeping bias voltage, are reported 
in Fig[31 in which only six of them, HOMO-2, • • • , 
LUM0-f2, are examined, since they dominate most of 
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the properties of the OSE during bias sweeping. The 
most striking feature of Fig[2] is the bias-induced en- 
ergy levels crossover. It occurs twice in the bias sweep- 
up [FiglHfa)] but only once in the sweep-down direction 
[FiglHb)]- The first crossover in FiglSja) does not in- 
crease the current. Detailed analysis finds it due to the 
fact that the corresponding HOMO and LUMO states 
are localized near the metal-OSE interfaces by the sur- 
face vibration modes (not shown here due to space limi- 
tation). Note that the HOMO/LUMO here is a reduced 
state in the presence of electrodes and lattice. Unlike 
the first crossover, the second one, shown by the dual- 
crossing between the HOMO and LUMO-l-1 and between 
the LUMO and HOMO-1 as indicated by the arrows in 
FigEl^a) around Vth, does increase the current, as it af- 
fects the system remarkably. It combines the valence and 
conductance bands of the OSE to form a single metallic 
band structure without energy gap; and meanwhile it also 
completely delocalizes the wave functions of those six lev- 
els over the whole OSE substructure The gap-free metal- 
lic band structure and the delocalized wave functions, 
together with the uniform valence bonds [cf. Figl^Jb)] 
around the second crossover, are highly advantageous to 
the transport of 7r-electrons. They lead to the I-M tran- 
sition of the organic electronics with rapidly increasing 
current around the Vth- 

In a similar but reversible manner, when the bias 
sweeps down the energy levels crossover at V ^ 2.8 V 
re-opens the band gap [FiglSl^b)], re- localizes the elec- 
tronic wave functions, and re-dimerizes the lattice dis- 
placements of the OSE. As results, the OSE substructure 
transforms from the metallic back to the insulating state 



(M-I transition) with the current sudden drop following 
a near-linear decrease (cf. FigdJ. 

Some experimental evidences of the bias-induced I-M 
transition have been reported in large size organic elec- 
tronics. One is the bias-driven high-to-low resistive state 
transition in the alkali-7,7,8,8-tetracyanoquinodimcthane 
(K-TCNQ) single crystals, a quasi-one dimensional or- 
ganic charge-transfer complex (Mott-Peierls insulator) 
The I-M transition occurs in this system at about 
several hundred of volt in low temperature, where a 
metallic path has been visualized with a microscope. The 
present study may also shed some light on understanding 
the experimental observations in Ref. |T^ . These include 
the conductance switching, hysteretic loop, and 'dimer- 
ization' reduction in the TTF-based dimeric donor salt; 
cf. Figll] here versus the Fig. 2 of Ref. [H. The detailed 
comparison between the theory and experiments will be 
published elsewhere. 

In summary, we have investigated the bias-induced I- 
M transition in organic electronics devices and explained 
it with the energy levels crossover which diminishes the 
Peierls phase and delocalizes the electron states at the 
threshold voltage. The bias-triggered intrinsic large on- 
off ratio in this work (> 20:1) is experimental control- 
lable. Our theory study reproduces a range of exper- 
imental observations, such as the TM phase transition 
(Til IT^ , the hysteretic conductance switching [l^, [S^l . 
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FIG. 1: The hysteretic current as the function of sweeping 
bias of the model metal/OSE/metal electronics at T = 11 K. 
The size of OSE is No = 40 and the linear voltage sweep rate 
is at 0.1 V/sec. The insert is the threshold voltage Vth as a 
function of the OSE chain length No- 
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FIG. 2: (color online) High (solid lines) and low (dash-dot 
lines) conductance states triggered by the up-sweeping bias, 
(a) Transmission coefficient T{E). The vertical dot lines in- 
dicate the window between and /ir for the charge trans- 
mission; (b) Lattice distortion ?/„ = — u„ of the OSE 
substructure. 
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FIG. 3: (color online) The voltage-evolution of six energy 
levels of the OSE (HOMO-2, • • • , LUMO+2) for up-sweeping 
(a) and down-sweeping bias (b). The arrows indicate the 
crossover of energy levels. 



